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Abstract 

Based on the world DIS data we extract the experimental F2 — F2 
as a function of Bjorken-x and photon virtuality Q^, using two differ- 
ent methods. Both methods lead to identical results. We find that the 
standard PDFs fail to describe the experimental data below < 7 
GeV^, which is much higher than for F2 and separately. The differ- 
ence between PDFs and the experimental data cannot be understood 
as due to nuclear effects in the deuteron, and evidently suggests sub- 
stantial nonsinglet higher-twist effects. The trend of the experimental 
data is approximately explained by a recent two-component model of 
the nucleon structure functions and suggests strong Q^-dependence of 
the Gottfried Sum, in disagreement with the parton model interpreta- 
tion. The large negative higher-twist effects can explain the difference 
between the value of the Gottfied integral obtained recently by the 
E866 Drell-Yan experiment at Fermilab and an older NMC result. 



The QCD-improved parton model (IPM) applies at large virtuality 
of the probe. However, the exact range of its applicability is not obvious. 
Generally a good desciption of i^f and F2 within IPM can be obtained 
already above ~ 1-2 GeV^ |, |, |. 

In the last decade the precision of experimental data on structure func- 
tions became so good that one can try to obtain the information on effects 
beyond the leading twists @, ^. 
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A recent lattice QCD calculation has found an unexpectedly large twist-4 
contribution to the structure functions |^] . This is rather difficult to reconcile 
in the light of the success of the improved parton model which is known to 
work empirically down to ~ 1-2 GeV^. 

In a recent work|^ we have shown that both the proton and deuteron 
structure functions can be well explained in a broad range of x and also 
in a two-component model of the structure function 

F^'^'ix, g2) = F2^/"'^"'^(:c, Q2) + Q2) , (1) 

which gives better agreement with data below ^ 2> GeV^. Our model 
fulfills F2^"'{x,Q'^) — > when — > 0, as required by current conservation. 
While ^"■''^"'^(x, Q^) ^ by construction, the vanishing of F2^"''^^^^{x, Q^) 
is achieved as in the low-x Badelek-Kwiecihski model B 



2^ 



where pP/^'°'^y"^P [g the standard parton-model structure function and x, 
are defined as in j^. This simple phenomenological form insures a correct 
Q"^ — >■ limit and is justified by the dispersion method For not too 
small X, i.e. in the region we are interested in here, the modification of the 
structure function arguments x —>■ x, Q'^ ^ Q'^ is not needed 

Except for a better description of the F2 data in the low-Q^ region the 
model from ||^] has some interesting consequences. The higher- twist contri- 
bution to the structure function can be defined by 

In the model for F2, F2 proposed in 0] there is a significant cancellation 
between positive (VDM-type) and negative (due to a modification of the 
partonic component) higher twist contributions 

5^^F2^/"(x,Q2) = 



Q + Qq 



Q + Qq 
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in the range of small and intermediate x. 

Neglecting a small p-uj nondiagonal term due to the exchange of a Regge 
trajectory 02 

we have 

where Cy is a normalization factor which can be related to the electro- 
magnetic V e^e~ decay width and Q{x,Q'^) is a form factor which phe- 
nomenologically accounts for the finite lifetime of the hadronic fiuctuation 
of the photon [|^. Therefore 

Dg{x,Q^) ^ Fi{x,Q')-F^{x,Q') = 



for the range of Bjorken-x we are interested in. Unlike the case of the 
structure functions there is no cancellation of hadronic and partonic higher 
twists. A strong Q^-dependence of Dq is obvious. The factor in front 
of the r.h.s. of Eq.(|^ assures also the vanishing of the Gottfried integral 
for Q2 _^ 0. This result is in contrast with the result obtained by Ball 
and Forte ^ who found a restoration of the Gottfried Sum Rule Sq = 
f {FKx, Q2) _ F^{x, Q2)} ^ 1 when ^ 0. 

In order to verify the anticipated effect discussed above let us come to 
the status of experimental data on F2 — F2 ■ Because F2 is not a directly 
measurable quantity one is forced to use rather deuteron data. The sim- 
plest approximation, neglecting all nuclear effects in the deuteron, would 
be to replace Dq = F^ - F^ by - {2Fi - F^). This direct method is, 
however, not very useful as far as the existing data are concerned. There 
the errors (both systematical and statistical) for F2 and F2 are indepen- 
dent and the resulting errors for Dq are too large to reveal the anticipated 
Q^-dependence. Instead of the direct subtraction of F2 and F2 one can use 
a method proposed by NMC at CERN In the following two methods 
will be used to extract Dq'- 
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. (A): (x, Q2) _ F^{x, Q^) = F^{x, Q^) ■ (l - , 

. (B): Ff (x,Q2) _ F2-(x,Q2) = 2Fi{x,Q^) ■ ^fg|} , 

where F2/F2 is extracted from the measured F2/F2 data. The latter was 
already used by NMC in the evaluation of the Gottfried integral in an un- 
published analysis of the dependence of the NMC data The two 
methods above are of course equivalent when nuclear effects are neglected. 
They are equivalent also in a less obvious case when nuclear effects are taken 
into account and F2/F2 is replaced by 2F2/F2 — 1 (this will be assumed 
hereafter). It can be shown that then both methods include nuclear effects 
Q in exactly the same way, identical as the difference - {2F^ - Ff). 
Both methods require knowledge of the ratio F2 / Fj for the same x and 
points as the structure functions F2 or F2 which is not the case for 
existing experimental data. For this reason use of a parametrization for 
F2/F2 appears necessary. In our analysis we have used a parametrization 
from Q. In Fig.l (top panel) we compare this parametrization with precise 
F2/F2 data from NMC The dashed lines indicate a 2 % uncertainty. 
Below we shall quantify the quality of the parametrization and justify the 
choice of its uncertainty band. In the bottom-left panel we present a per- 
centage of points above (solid line) and below (dashed line) a trial (rescaled) 
parametrization for the ratio, as a function of the relative deviation from 
the original parametrization. The result clearly indicates that statistically 
the same amount of points lies above and below the nominal curve. In the 
bottom-right panel we show the fraction of NMC data points in a band of 
uncertainty as a function of the band width (dashed line). About 3/4 of 
the experimental points remain in the 2 % band of uncertainty. Because the 
functional form of the ratio is not known, and some fluctuations of the ratio 
are not excluded a priori, we have also performed the following analysis. 
We assign to each experimental point a Gaussian distribution of probability, 
centered on the experimental point, with a standard deviation equal to the 
experimental error. We then determine the fraction of the area that lies be- 
low such a probability distribution within the parametrization uncertainty. 
We show this quantity averaged over 260 NMC experimental points [^] in 
the bottom-right panel by the solid line marked by < P >. About 60 % of 

^the latter will be marked by a star in Fig. 2 and 3. 
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so-defined probability remain in the 2 % band of uncertainty. 0. This anal- 
ysis assumes inherently independence of successive experimental points and 
allows fluctuations around the smooth parametrization. Because we do not 
expect neither a full independence of neighbouring experimental points nor 
sharp fluctuations of the ratio, the 2 % uncertainty band seems to be rather 
an upper limit for the uncertainty on F2/F2. Taking into account uncer- 
tainty of F^/F^ at its upper limit should allow us to prove our thesis about 
the IPM breaking with a bigger degree of reliability. We have checked that 
in the interesting kinematical range the used parametrization is consistent 
within experimental uncertainties with the E665 [^] data for F2 /F2. 

The world data on F2 and F2 from the compilation was taken in 
the present analysis. 

In Fig. 2 we show experimental data for Dg{x) obtained with the two 
methods for three different values of photon virtuality = 1.1, 3.4, 7.0 
GeV^ as a function of Bjorken-x. The error bars take into account the 
inherent uncertainties of -Ff (method A) or F2 (method B) as well as the 
uncertainty of F2/F2. The two bands below each distribution separately 
show the errors due to the mentioned above 2% uncertainty in F2 /F2 and 
due to the uncertainty of F2 (Ihs) or F2 (rhs). As seen from the figure the 
error of Dq due to the uncertainty of the F2 / F2 ratio is substantially bigger 
than the one due to the uncertainty of either F2 or i<2^. Within experimental 
errors the two methods lead to identical results despite the different data sets 
used. We have also checked that the result obtained by a direct subtraction 
of the deuteron and proton structure functions leads to a consistent result 
although the scatter of the data is much larger. Therefore the data obtained 
in this way can be reliably used to study the dependence of F2 — F2 
and/or corresponding higher- twist effects. 

For comparison we show theoretical results obtained with GRV94 p] 
(short-dashed line), MRST ||2| (long-dashed line) and a recent CTEQ5 ||] 
(solid line) NLO DIS-scheme parton distributions. ^ Above > 7 GeV^ the 
parton model describes the experimental data well. However, the predictions 
of the improved parton model clearly deviate from the experimental data 
below <7 GeV^ for Bjorken x between 0.2 and 0.45. This deviation is not 
a random fluctuation. It is present in a broad range of as will be shown 
in the next figure. Similar deviations of the IPM from the experimental data 

^Please note that within standard statistical analysis, assuming the Gaussian proba- 
bility distribution, the required fraction for the standard deviation is 68.2 % 

^The MRST result is shown only above > 1.25 GeV^ and CTEQ5 above > 1 
GeV^, i.e. in their applicability range. 
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can also be observed for F2 and F2, albeit at smaller values of This 
is consistent with the cancellation of higher-twists for F2 or F2 and with the 
lack of such a cancellation for F2 — F2 as discussed above. None of nuclear 
effects we know can explain this disagreement. Moreover the agreement 
of the improved parton model with experimental data for small values of x 
is most probably accidental, as shadowing effects which are not included in 
this calculation, but are present in the data, would enhance the theoretical 
results by 0.005-0.015, depending on (x, Q"^) and the model used. This would 
worsen the agreement of the IPM also for small values of Bjorken x < 0.1. 
For the sake of simplicity, we omit nuclear corrections in the present paper. 
We have estimated them and found that they do not affect our conclusions. 
For comparison we show the prediction of our model Q (thick solid line) 
with Ql = 0.79 GeV^, with target mass corrections calculated according to 
[16|. The target mass effects become important at small and larger x. 

In order to better visualize the unexpected within the IPM Q^-dependence 
we present in Fig. 3 Dq as a funtion of for some selected values of Bjorken 
X = 0.140, 0.275, 0.350. This figure clearly demonstrates a failure of the 
IPM as well as substantial negative higher- twists for F2 — -Pf . The trend of 
the experimental data is completely different to the prediction of the stan- 
dard IPM. Our model reproduces well the main trend of the experimental 
data and provides a better description than the standard IPM. However, 
the agreement of our simple model with experimental data is not complete. 
The structure function difference F2 — F2 is a relatively small quantity and 
therefore may be sensitive to subtle effects. At low , i.e. rather low energy 
for fixed Bjorken-x, one should worry about a careful treatment of different 
exclusive channels and their contribution to F2. One may also expect some 
contribution due to the p — uj nondiagonal VDM term NNLO analysis 
and careful treatment of nuclear effects could also be valuable. 

The effect discussed in the present paper may be very important to un- 
derstand small differences for light-antiquark flavour asymmetry obtained 
from different types of experiments. Recently the E866 experiment on Drell- 
Yan production in proton-proton and proton-deuteron scaterring has pro- 
vided a new, rather precise, information about the nucleon d — u asymmetry 
[p^ . Because of the broad range of x available in this experiment we could 
estimate ^ 

I dx [d{x) - u{x)]Esm = 0.09 ± 0.02 . (8) 
h 

No strong Q^-dependence is predicted for nuclear effects except of a VDM-type 
shadowing for x < 0.1. 
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This number is smaller than the one deduced from a somewhat older NMC- 
result lie] on the Gottfried Sum RuleQ 



dx [d{x) - u{x)]nmc = 0.147 ± 0.039 . (9) 

JO 

While the average photon virtuality in the Fermilab experiment was fairly 
large < >^ 50 QeV^, the average photon virtuality in the CERN deep 
inelastic scattering experiment was relatively low: < >^ 4 GeV^. The 
latter experiment is in the range of where we have observed large negative 
higher-twist effects (see Figs. 2 and 3). This means that the deviation of the 
Gottfried Sum Rule from its classical value of ^ ||l^ observed by NMC [^] 
is partially due to the d — u asymmetry and partially due to the discussed 
higher- twist effects as follows from Eq.(|^. The higher- twist effects would 
resolve then the discrepancy between (^) and (P). 

In summary, we have found substantial deviations from the improved 
parton model for the difference F2 — F2 already at as large as 7 GeV'^, 
which is much higher than for F2 and F2 separately. The deviation has 
been predicted by us in Q. When combined, these two analyses strongly 
suggest that the agreement of the improved parton model for F2 and F2 at 

< 7 GeV^ is rather accidental. The present result is in contrast to rather 
small higher- twist effects estimated in F^^ {x,Q'^) The nuclear effects 
in the deuteron have not been included in the present analysis. However, 
we estimated their role within simple nuclear models and found that tak- 
ing into account these effects would only strengthen our conclusions. The 
results obtained here strongly suggest the Q^-dependence of the Gottfried 
Sum Rule and could resolve a discrepancy between the u — d asymmetry 
obtained in different processes. The breaking of the parton model approach 
for the reported range would have important consequences for different 
analyses of exclusive reactions concerning the flavour and spin structure of 
the nucleon where the applicability of the parton model is conditio sine qua 
non. 



^We note that both analyses assume isospin symmetry for proton and neutron quark 
distributions. 
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Figure 1: The quality of the F2/F2 ratio parametrization for the NMC 
F2/F2 data ^J^/. In the top panel shown is F2/F2 shifted by 0, 0.15, 
1.5 for X = 0.675, 0.4-50, 0.0015, respectively, as a function of the pho- 
ton virtuality. The short-dashed line shows a 2 % uncertainty of the ratio as 
discussed in the text. In the bottom-left panel shown is an up- down asym- 
metry of the NMC data with respect to the used parametrization and in the 
bottom-right panel a fraction of the same data points in the parametrization 
uncertainty band. 
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Figure 2: Dq{x) for some selected values obtained with the method A 
(Ihs) and the method B (rhs). Below each distribution shown separately is an 
error band due to the uncertainty of the F2 jF^ ratio and an error band due 
to the uncertainty of Fl^ (Ihs) or F2 (rhs). The short-dashed line correponds 
to the GRV94 NLO PDF's, the long-dashed line to the MRST98 NLO 
^ PDF's and the solid to the CTEQ5 NLO Q PDF's. The thick solid line 
represents a prediction of our model Jl/ with details as described in the text. 
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Figure 3: Dg{Q^) for some selected values of x. The meaning of the curves 
is the same as in Fig. 2. 
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